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Abstract. The acceleration parameter defined through the local volume expansion is 
negative for a prcssureless, irrotational fluid with positive energy density. In the presence 
of inhomogeneities or anisotropies the volume expansion rate results from averaging over 
various directions. On the other hand, the observation of light from a certain source 
in the sky provides information on the expansion along the direction to that source. If 
there are preferred directions in the underlying geometry one can define several expansion 
parameters. We provide such definitions for the case of the Tolman-Bondi metric. We 
then examine the effect of a localized inhomogeneity on the surrounding cosmological 
fluid. Our framework is similar in spirit to the model of spherical collapse. For an 
observer in the vicinity of a central overdensity, the perceived local evolution is consistent 
with acceleration in the direction towards the center of the overdensity, and deceleration 
perpendicularly to it. A negative mass leads to deceleration along the radial direction, 
and acceleration perpendicularly to it. If the observer is located at the center of an 
overdensity the null geodesies are radial. The form of the luminosity distance as a 
function of the redshift is consistent with acceleration for a certain range of redshifts. 
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1. Introduction 

The cosmological expansion of our Universe seems to have accelerated in the recent past. 
This conclusion is supported by the form of the luminosity distance as a function of the 
redhshift for distant supernovae [1, 2]. In the context of homogeneous cosmology, a recent 
accelerating phase is also in agreement with the observed perturbations in the cosmic 
microwave background [3]. The mechanism that triggered the acceleration has not been 
identified conclusively. The simplest explanation is that the cosmological constant is non- 
zero. However, the absence of acceleration at redshifts z ^ 1 implies that the required 
value of the cosmological constant is approximately 120 orders of magnitude smaller than 
its natural value in terms of the Planck scale. 

One intriguing fact is that the accelerating phase coincides with the period in 
which inhomogeneities in the matter distribution at length scales & 10 Mpc become 
significant, so that the Universe cannot be approximated as homogeneous any more at 
these scales. A link between inhomogeneities and cosmological acceleration has been 
pursued in various studies. There have been arguments, based on perturbative estimates, 
that the backreaction of superhorizon inhomogeneities on the cosmological expansion is 
significant and could cause the acceleration [4]. However, the validity of this effect is 
questionable [5]. 

We are interested in the importance for the problem of cosmological acceleration of 
inhomogeneities with sub-horizon characteristic scales today. Because of the significant 
growth of such inhomogeneities at recent times, a perturbative treatment may not be 
sufficient. An exact solution of the Einstein equations, even for a simplified geometry, 
could be more useful in order to reveal an underlying mechanism. The Tolman-Bondi 
metric [6] has been employed often in this context [' -[12]. It has been observed that any 
form of the luminosity distance as a function of the redshift can be reproduced with this 
metric [10]. 

A drawback of the standard handling of the Tolman-Bondi metric is that the radial 
coordinate is defined such that the fluid density is initially perceived as homogeneous. 
The presence of inhomogeneities is introduced through a function that determines the 
local Big Bang time. This obscures the intuition on the role of large mass concentrations. 
Morever, the growth of perturbations and its effect on the expansion is not obvious. 

In our analysis we choose a gauge such that the initial density perturbation is 
apparent. We first model the perturbation by matching a Schwarzschild metric in 
the interior with an exterior Tolman-Bondi metric. We show that the presence of 
a large overdensity modifies the cosmological evolution of the surrounding fluid. In 
particular, accelerating expansion can take place along the radial direction and can be 
observed through the redshift of light signals propagating radially. On the other hand, 
the expansion perpendicularly to the radial direction is characterized by deceleration. 
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Underdensities do not induce acceleration in any direction. 

As an exotic possibility, we also consider the case of a negative mass at the center. 
A fluid with energy density that can become negative has already been discussed in 
the literature. In ref. [13] a field termed phantom with negative kinetic energy drives 
the accelerating expansion. Our setup is based on less drastic assumptions. We consider 
vacuum solutions of the Einstein equations that can be interpreted as arising from objects 
of negative mass. The relevant metric is the Schwarzschild one with a negative mass. We 
show that a pressureless fluid near the center of such a configuration can have accelerating 
expansion perpendicularly to the radial direction. 

We should mention at this point that negative-mass configurations have several 
features that appear problematic, and their physical significance is not guaranteed. The 
negative-mass Schwarzschild solution has a naked singularity at the center, so that it 
violates the cosmic censorship conjecture. However, naked singularities appear often in 
studies of gravitational collapse for large classes of initial conditions [I I]. Another issue 
concerns the stability of the negative-mass Schwarzschild solution, which has not been 
proven in full generality. In ref. [16] the linear stability has been demonstrated, under a 
family of boundary conditions at the singularity that have a physical motivation. On the 
other hand, runaway solutions of the Einstein equations are known for a pair of objects 
with positive and negative mass [15, 16]. We do not offer here a resolution of these 
problematic issues. Instead, we point out an interesting feature of the negative-mass 
Schwarzschild solution related to the issue of the cosmological acceleration. 

Next we consider a homogeneous energy distribution in the central region. We define 
an initial condition for the spatial curvature which makes our model very similar to that 
of spherical collapse [17]. A small overdensity grows consistently with the perturbative 
Jeans analysis, until the density fluctuation becomes comparable to the average density. 
Then, the overdensity starts collapsing. We study the null geodesies that start outside the 
perturbation and lead to an observer in its interior. If the observer is located in the center 
of the spherical perturbation the geodesies follow the radial direction. The accelerated 
expansion of the surrounding fluid is reflected in their form, which, in turn, determines 
the luminosity distance of a light source as a function of the redshift. 

We demonstrate that the expected luminosity distance is consistent with accelerated 
expansion for a certain range of redshifts. In this work we do not put emphasis on the 
exact numerical consistency with the data. The reason is that we model the fluid in 
the central collapsing region as pressureless. This results in a strong deceleration for 
small redshifts, which is not expected to be present if the cosmological fluid can reach 
virialization. Our main point is the identification of a mechanism that could drive the 
observed acceleration. This mechanism is based on the gravitational force that acts on 
the fluid surrounding a collapsing overdensity and its effect on the expansion rate. 

In the following section we discuss the cosmological expansion in inhomogeneous 
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cosmologies and its dependence on the direction of observation. In section 3 we employ 
the Tolman-Bondi metric in order to construct a model for a central inhomogeneity in 
an asymptotically homogeneous Universe. In section 4 we examine if acceleration can be 
induced by the inhomogeneity along or perpendicularly to the radial direction. In section 
5 we study the form of the luminosity distance as a function of the redshift for an observer 
located at the center of an overdensity. In section 6 we provide a summary of the results 
and our conclusions. 



2. Acceleration in inhomogeneous cosmology 

The presence of a center in the configuration that we would like to study implies that we 
have to assume a metric appropriate for an inhomogeneous fluid. Under the assumption 
of spherical symmetry, the most general metric for a pressureless, inhomogeneous fluid is 
the Tolman-Bondi metric [()]. It can be written in the form 

ds 2 = -dt 2 + b 2 (t, r)dr 2 + R 2 (t, r)dn 2 , (2.1) 

where dQ 2 is the metric on a two-sphere. The function b(r, t) is given by 

where the subscript denotes differentiation with respect to r, and f(r) is an arbitrary 
function. The bulk energy momentum tensor has the form 

T A B = diag(-p(t,r), 0, 0, 0) . (2.3) 

The fluid consists of successive shells marked by r, whose local density p is time-dependent. 
The function R(t, r) describes the location of the shell marked by r at the time t. Through 
an appropriate rescaling it can be chosen to satisfy [14] 

R(0,r)=r. (2.4) 

The Einstein equations reduce to 

M, r {r) = 4irR 2 pR >r , (2.6) 

with G = (167rM 2 ) . The generalized mass function M.{r) of the fluid can be chosen 
arbitrarily. It incorporates the contributions of all shells up to r. It determines the energy 
density through eq. (2.6). Because of energy conservation M.{r) is independent of t, while 
p and R depend on both t and r. 

The volume expansion rate is usually defined through the four-velocity of the fluid 
u a as 

3H = u\ = Ua , b g ab = u a , b h a \ (2.7) 
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where 

h ab = g ab + u a u b_ ^g) 

The Raychaudhuri equation constrains its evolution with time. It reads 

H = H >c u c = -H 2 - l -a ab a ah + ^ ab cu ab - ^R ab u a u\ (2.9) 

where a ab and uo ab are the shear and vorticity tensors respectively. For an observer 
comoving with a pressureless and irrotational fluid, the acceleration parameter is 

" = ip + 1 = -h (r° tff0i + mp f ) ■ (2 - 10) 

If the local energy density is positive, no acceleration can take place. 

In the context of inhomogeneous cosmology the definition (2.7) does not capture the 
variation of the expansion rate in different directions. For example, the Tolman-Bondi 
(2.1) has a preferred direction (the radial one). In such situations we can define a tensor 
p ab that projects every quantity perpendicularly to the preferred space-like direction s a 
(and of course the time-like vector field u a ). This is 

f h = g ab + u a u b - s a s b = h ab - s a s b . (2.11) 

For the Tolman-Bondi metric (2.1) we have s a = b~ l d r . We can now define invariant 
expansion rates parallel and perpendicularly to s a , according to [8] 

H r = u a ., b s a s b = b -± = ^ (2.12) 

H e = l -u a , b p ab = ^ ) (2.13) 

so that 

H= 2 -H e + l -H r . (2.14) 

From the above, it is obvious that H corresponds to an average of the expansion rates 
in various directions. If the acceleration is defined through H, the expansion is always 
decelerating for p > 0. This is the conclusion drawn in previous studies [7, 8, 9]. However, 
it is possible that the expansion may be accelerating is some direction even though eq. 
(2.10) always holds. 

A definition of the expansion rate that takes into account its directional dependence 
is given in ref. [19]. It is 

H= l -u a . a + a ab J a J b , (2.15) 

where a ab is the shear tensor and J a a unit vector pointing in the direction of observation. 
In the context of the Tolman-Bondi metric, and for an observer located away from the 
center of the configuration, the above definition gives [19] 
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The parameter tjj is the angle between the radial direction through the observer and the 
direction of observation. For ip = or 7r we have H = H r , while for if; = n/2 or 3ir/2 we 
have H = Hq. 

A definition of the acceleration parameter in a specific direction can be given in terms 
of the expansion of the luminosity distance Dl of a light source in powers of the redshift 
z of the incoming photons. For small z this is 

J 2 D, 
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dz 2 
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In ref. [19] the acceleration parameter was calculated in the context of the Tolman-Bondi 
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These parameters do not depend solely on local quantities, as opposed to the acceleration 
parameter of eq. (2.10). For example, from eq. (2.5) we have 

1 Ml , 

qe = -Tmp&h 2 - (2 - 25) 

The value of qg depends on the total mass function and not just the local energy density. 
For M > we expect to observe deceleration perpendicularly to the radial direction, 
while for M < acceleration. 

It must be emphasized that the expressions (2.23), (2.24) refer to the effective 
acceleration deduced from light signals that originate in the vicinity of an off-center 
observer. In this sense, they determine the character of the local expansion. It is 
interesting that they do not predict a dipole component [19]. On the other hand, light 
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signals that travel longer distances (e.g. from outside a spherical inhomogeneity to its 
center) integrate information on the local expansion along their whole trajectory. We 
shall discuss the form of the luminosity distance as a function of the redshift for this case 
in a following section. 

Finally, we point out that in the standard treatment of the Tolman-Bondi metric 
-[12] the radial coordinate is taken such that the initial energy density is perceived 
as homogeneous. Moreover, the function f(r) is assumed to be zero in order to provide 
consistency with a matter dominated, flat Universe. The inhomogeneity is introduced 
through a function tb(r) that appears in the integration of the Friedmann equation and 
determines the local Big Band time. Our approach is closer to the physical picture, as the 
overdensities are apparent from the beginning. Moreover, we put the emphasis on their 
effect on the cosmological evolution of the surrounding fluid. 



3. Expansion near an inhomogeneity 

We would like to describe the cosmological evolution around a central region in which the 
mass function deviates from the form implied by a homogeneous distribution of matter. 
We have in mind the general form 

M(r) = M + SM(r), (3.1) 

with M.q positive or negative. Even for negative M.q we assume that eq. (2.6) is satisfied 
with p > 0. No negative energy density appears in the energy- momentum tensor, and the 
weak energy condition holds. We now have 

8M, r =AixS 2 pS )T . (3.2) 

The most convenient way to realize our scenario is by matching two spherically 
symmetric space-times. We assume that in the central region the metric takes the 
Schwarzschild form with mass Aio 

ds 2 = -A 2 (f)dP + B 2 (f)df 2 + fW, (3.3) 

where A 2 (r) = B~ 2 (f) = 1 — M.q/ {&irM 2 f). The Tolman-Bondi metric (2.1) determines 
the geometry at large distances. The two regions are separated by a spherical surface 
with metric 

ds 2 = -dr 2 + S 2 (j)dVL 2 . (3.4) 

The physical situation described by our scenario does not involve necessarily a singularity 
at the origin. For Ai > the metric of eq. (3.3) may result from a compact object 
with positive energy density and radius smaller than S(t). For Ai < we are forced to 
consider eq. (3.3) as a vacuum solution in order to avoid a negative energy density. 
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The matching of (2.1) and (3.4) leads to 

S 2 (r) = R 2 (t(r),r(r)) (3.5) 

-i 2 (r) + b 2 (t(T),r(r)) r 2 (r) = -1, (3.6) 

where a dot denotes a derivative with respect to the proper time r on the surface. The 
matching of (3.3) and (3.4) gives 

S 2 (r) =f 2 (r) (3.7) 

-A 2 (f(r)) t (r) + £? 2 (f(r)) rV) = - 1. (3.8) 

The absence of shell crossing implies r(r) = 0. The shell that coincides with the matching 
surface is always the same. As a result, it maintains a constant value r(r) = r . Eq. (3.6) 
then gives £(r) = r, while eq. (3.5) can be written as S(r) = R(r,r ). 

In order to obtain an equation for S(t) we need the additional assumption that there 
is no singular energy density on the matching surface. This implies that the extrinsic 
curvature is continuous across the surface. The resulting condition reads 

s{T) [ s{T) + WsW)J = — W^) — ' (3 ' 9) 

Making use of eq. (2.2) we find 

*M = '™ + ssw (X10) 

This coincides with the equation of motion of the first shell, as derived from eqs. (2.5), 
(2.6). 

We can conclude that, in the absence of shell crossing, the interior Schwarzschild 
metric can be matched with the exterior Tolman-Bondi one on a spherical surface that 
coincides with the location of the first shell S(t, r ). The equation of motion for the shells 
of fluid is 

R 2 (r,r) 1 M(r) f(r) 

R? 8ttM 2 R 3 R 2 ' 1 ' ' 

where A4(r) is given by eq. (3.1) with 8M,(r) — > for r — ► r$. The initial condition 
is given by eq. (2.4). The above equation can be considered as a generalization of the 
standard Friedmann equation for an inhomogeneous cosmological fluid. The integrated 
mass within a spherical volume of comoving radius r is given by Ai(r), while the function 
f(r) defines a generalized curvature term. 

In order for eq. (2.5) to have a solution, the two arbitrary functions f(r) and «M(r) 
must satisfy f(r) > —Ai(r)/(87rM 2 r) for all r. The function f(r) defines an effective 
curvature term in eq. (2.5). We can also interpret f(r) as part of the initial radial 
velocity of the fluid. This has to be non-zero in the presence of a negative mass for a 
physically meaningful solution to exist. 
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As a physically motivated example we consider a fluid with energy density that is 
initially constant for r > ro, so that 5M.(r) = Anpo (r 3 — r$) /3. We assume that a 
spherically symmetric inhomogeneity has appeared during the cosmological evolution of 
the fluid (probably as a result of gravitational instability). At tj = the inhomogeneity 
is concentrated in the region r < r$. We parametrize the total mass of the inhomogeneity 
as Ai = 47rprQ/3 with p = ep . For e > 1 the inhomogeneity is a local overdensity, for 
< e < 1 an underdensity, while for e < we have an object of negative mass at the 
center of the configuration. 

Our crucial assumption is that the expansion rcitc £it some initial time U is given for all 
r by the standard expression in homogeneous cosmology: Hf = po/(6M 2 ). We essentially 
assume that the expansion is completely homogeneous at the time of appearance of the 
inhomogeneity. This is very similar to the initial condition considered within the context 
of the model of spherical collapse [17, 18]. Then, eq. (3.11) with R(0,r) = r implies that 
f{f) = (1 — e)/f, where / = f/Hf and Hi = Hir . An effective curvature term appears 
that is proportional to the size of the inhomogeneity e — 1. For f — > oo we have / — > 0, so 
that the curvature term disappears and we recover the standard expansion in a spatially 
flat Universe. 

Eq. (3.11) can be written as 

-2 e — 1 + f 3 1 — e 

R(r,f)= + +— , (3.12) 

where f = tHi = tH{, R = R/tq, f = t/tq, and the dot denotes a derivative with respect 
to f. Sufficiently far from the center of the configuration (large f), eq. (3.12) reduces to 
the standard Friedmann equation for a homogeneous, matter dominated, flat Universe: 
(R/Rj ~ R~ 3 . Near the center of the configuration the expansion is inhomogeneous. 

For the value e = 1 we recover the standard homogeneous expansion. In this case the 
total mass in the region f < 1 could arise from a fluid with constant initial energy density 
equal to p$. The range e > 1 corresponds to an overdensity in the region f < 1. This 
could result from a continuous distribution of matter, or from a singularity at the center, 
such as a black hole. The range < e < 1 corresponds to an underdensity. Finally, the 
range e < corresponds to a negative mass at the center. 

In the case of an overdensity (e > 1), there is a certain time f c when R(f c , 1) is 
sufficiently large for the r.h.s. of eq. (3.12) to vanish at f = 1. At later times the first 
shell must reverse its motion and contract. This is very similar to the phase of collapse 
in the spherical collapse model. 



Cosmological Acceleration and Gravitational Collapse 



10 




r 

Figure 1. The solution R(f,f) of eq. (3.12) for e = 1.5. 

4. The local acceleration 

4-1. Positive mass 

In fig. 1 we plot the solution of eq. (3.12) for e = 1.5, with the initial condition R(0, f) = f. 
We observe deviations from the initial linear relation between R(f, f) and f. These are 
induced by the overdensity at f < 1. In particular, the additional gravitational attraction 
on the shells with f ~ 1, resulting from the overdensity, slows down the outward motion of 
these shells. At times f ~ 6 the first shell is beginning to reverse its motion and collapse 
towards the center. The fluid far from the overdensity is not affected significantly by it. 
Its expansion is typical of a homogeneous fluid. 

The evolution in fig. 1 is very similar to the one observed in the spherical collapse 
model [17]. An initial overdensity decouples from the homogeneous expansion and 
eventually reverses its motion and collapses. Our model provides additional information 
on the evolution of the material that surrounds such an overdensity. In particular, the 
relative motion of two fluid shells can be determined. It is apparent from fig. 1 that the 
f-derivative of the function R(f,f) is largest near f — 1. This implies that the relative 
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Figure 2. The local acceleration parameter q r {T,r) in the radial direction for e = 1.5. 



separation of the shells nearest to the center increases in comparison to the unperturbed 
fluid at large f. For a comoving observer with f ~ 1 the expansion seems faster than in 
the homogeneous model. On the other hand, two observers with the same radial distance 
from the center, but located at different angles, have a relative separation that tends to 
increase at a slower rate than for a homogeneous fluid. 

This behaviour is reflected in the values of the acceleration parameters. They can be 
expressed as 
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(4.2) 



where b = R' /J + f, f = f /Hf, Hi = Hiro, the dot denotes a derivative with respect 
to f , and the prime a derivative with respect to f. We use Hi = 1 throughout this 
section. The parameter q r (f,f) in the radial direction is plotted in fig. 2. Initially it is 
negative for all f. At later times it becomes positive in the central region. At all times 
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Figure 3. The expansion rate H and the local acceleration parameter q, as a function 
of the angle ip relative to the radial direction, for e = 1.01, f = 1 and f = 100. 



and for large f it approaches asymptotically the value q r = —0.5, typical of a homogenous 
matter-dominated flat Universe. The acceleration parameter qg(f,f) perpendicularly to 
the radial direction is always negative. For large f it approaches the value q$ = —0.5 as 
well. 

The appearance of accelerating expansion does not require e>l necessarily. In fig. 
3 we plot the expansion rate H and the acceleration q as a function of the angle ip for a 
model with e = 1.01. For this we employ the expressions (2.16) and (2.18)-(2.22). We 
consider the first shell with f = 1 at a time f = 100. The dependence of the expansion 
rate on the angle is mild. We can observe signs of the slowing down of the relative 
expansion perpendicularly to the radial direction: Hq = H(n/2) < H r = H(Q). On the 
other hand, the system is far from the period of collapse, which will take place at times 
f = (9(10 3 ). At f — 100 the ratio of the average density in the region r < ro and the 
density at large r is po/Poo — [(-Roo/ r oo)/(-Ro/ r o)] 3 e — 1.18. The local overdensity can still 
be considered a perturbation of the homogeneous background. Despite that, it induces 
an effective acceleration close to the radial direction, for angles ip ^ tt/6. 

For < e < 1 both acceleration parameters remain negative at all times. This means 
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that the presence of a void cannot induce accelerating expansion in the neighbouring 
regions. 

4-2. Negative mass 

In fig. 4 we plot the solution of eq. (3.12) for e = —0.5, with the initial condition 
R(0,f) = f. As in fig. 1, we observe deviations from the initial linear relation between 
R(f, f) and f. In this case, however, the shells nearest to the center are repelled strongly 
by the negative mass. They expand faster than the unperturbed medium at large f. At 
times f ~ 5 the first shell catches up with the adjacent one (the f-derivative of R(f,f) 
vanishes at f = 1) and the phenomenon of shell crossing appears. This is caused by the 
absence of pressure in our approximation. Our model is not adequate for the discussion 
of the evolution at later times. 

For < f ^ 1 the relative separation of adjacent shells near f = 1 increases, 
but with a decreasing rate. This is reflected in the negative value of the acceleration 
parameter q r (f,f). For f ^ 1 the relative separation decreases. As a consequence, an 
observer comoving with the fluid should perceive contraction in the radial direction. This 
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is confirmed by the value of the expansion rate H r of eq. (2.12), which becomes negative 
for f ^ 1 near f = 1. On the other hand, the relative separation of two observers with 
the same f, but located at different angles, increases faster near f = 1 than for large f. 
The acceleration parameter qo(f,r) perpendicularly to the radial direction is plotted in 
fig. 5. From eq. (3.12) we find 

al 2 R e-l+f 3 fA . 

= Si?-- (43) 
This implies that qe(f,r) is positive for f < 1.5 1 / 3 ~ 1.15 at all times for our choice of 
parameters. For large f both acceleration parameter approach asymptotically the value 
q = —0.5, typical of a homogenous matter-dominated flat Universe. 



5. The luminosity function for a central observer 

In this section we would like to study the issue of accelerated expansion for an observer 
within the perturbation. We assume that the observer is located at the center of a 
spherically symmetric overdensity and receives light signals from distances that extend 
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Figure 6. The null geodesic r g (f) and the redshift z(f) for a model with e = 1.01, 
Hi = 1, f e = 1461, f = 1563. 

beyond its surface. The assumption of spherical symmetry makes the problem tractable. 
In a realistic scenario we expect deviations from the exact symmetry, but the essense of 
the mechanism should remain unaffected. On the other hand, the preferred location of 
the observer is a more profound assumption. It could be argued that we are located near 
the center of a galaxy, which is a significant perturbation in the average density. However, 
it is questionable if an effect at the galactic scale can leave traces in the expansion at the 
level of the horizon. On the other hand, the mechanism we consider may be relevant 
for inhomogeneous cosmologies with many centers. The study of this situation is very 
difficult technically, and some kind of averaging will have to be implemented (see e.g. 
[20]). The necessary formalism in order to address this problem has not been developed 
yet. 

The other important feature of the scenario with an observer at the center of 
symmetry is that all the incoming signals follow radial geodesies. The radial acceleration 
can become positive, in contrast to the tangential one. In this sense, the central location 
of the observer provides the optimum scenario for our purposes. 

We also mention that our assumption about the location of the observer provides 
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automatic consistency with the isotropy of the cosmic microwave background. Deviations 
from complete isotropy, as quantified in the difference of power between the northern and 
southern hemisphere or possible allignments of the low multipoles [21], could possibly 
allow for an off-center location as well. An alternative possibility is that the isotropy 
appears at a sufficiently large length scale in a configuration with several centers of 
symmetry. It remains to be seen if the radial acceleration can be the dominant 
contribution for a random location of the observer. 

The simplest example of a perturbation on a homogeneous background has 

p = epo for r < r (5.1) 

p = Po for r > r . (5.2) 
The curvature term is 

/ = (l-e)f 2 forf<l (5.3) 

f=^—^- forf>l. (5.4) 

r 

For e > 1 the inhomogeneity is a local overdensity, for < e < 1 an underdensity, while 
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Figure 8. The acceleration parameter q and the local radial acceleration q r as a function 
of the redshift z, for a model with e = 1.01, Hi = 4, tq = 1567, zq = 0.067. 



for e = 1 we have a homogeneous energy distribution. The Friedmann equation becomes 

;3 ^2 



^2 

R r" r" 
R e-l + f 3 1 - e 



for r < 1 (5-5) 



# i?3 ' fiP f ° rf>1 ' (5 ' 6) 

where the dot denotes a partial derivative with respect to r = tHi. For e = 1 we recover 
the standard equation for a homogeneous, matter dominated, flat Universe. For e > 1 
the interior (f < 1) expands as a homogeneous, closed Universe with spatial curvature 
proportional to e — 1. The exterior (f > 1) is inhomogeneous. For large f we recover the 
Friedmann equation for a homogeneous, matter dominated, flat Universe: 

(5.7) 



-2 

R 



r 

R 2 W 

We can obtain an analytical expression for the growth of an overdensity in our model. 
We define the quantity 

c 



C(f) 



R(f, 1) 



1. 



(5.1 
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with ^> r = 1. The ratio of the energy density within the perturbation to the energy 



density far from it is given by the factor 
and f <: 1 

6 /3 x 2 ' 3 



C * 5 [f 



[1 + C) 3 - Using eqs. (5.5), (5.7) we find for ( ^ 1 



(5.9) 



where S = e — 1. At a time f 2 ~ 5~ 3 / 2 we have (1 + Q 3 ~ 2, which means that the 
energy density within the perturbation is double the asymptotic value. The phase of 
gravitational collapse starts at a time f c ~ 1.5 r 2 . The growth of a perturbation ~ f 2 / 3 is 
in qualitative agreement with the behaviour predicted by the Jeans analysis for subhorizon 
perturbations in a matter dominated Universe. It is also consistent with the growth of 
superhorizon perturbations in the matter dominated era. Our model can be viewed as an 
exact solution of the Einstein equations that is consistent with the behaviour expected 
from perturbation theory at all the scales of its applicability. This implies that we can 
consider values larger that one for the dimensionless quantity Hi = HiTo, for which the 
initial perturbation extends beyond the horizon. 

On radial null geodesies we have ds 2 = dVt 2 = 0, so that the geodesic equation 
becomes 

1/2 

(5.10) 



dr g (r) 
df 



H- 2 + f(f 9 (T)) 



where we have considered an incoming signal and the prime denotes a partial derivative 
with respect to f. We first determine numerically a solution of eqs. (5.5), (5.6) and then 
use the function R(f,f) in order to derive a solution of eq. (5.10). 
The redshift is given by the expression [10] 



lnfl 



R [f(f 



9)i 



[1 + / ft)] 



^-df 

1/2 "'S 



f0 R [f,f g { f )\ 

em R' [f,fg(f)} 



df. 



(5.11) 



The signal is emitted at a time r em from a point with comoving coordinate r em , and is 



received at the present time r at r 



D L = {1 + z) 2 R(f ( 



emi ' em 



0. The rescaled luminosity distance is [11, 19] 

(5.12) 



Through variation of r em we obtain the function Dl{z) = r^D^z). The expansion rate 
and the acceleration as a function of z can be defined as [1] 
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(5.13) 
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where the primes denote derivatives with respect to z. 

In fig. 6 we plot the null geodesic f g (f) for a model with e = 1.01 corresponding to 
a small initial overdensity at the level of 1%. The initial time is f, = 0. We take Hi = 1, 
so that initially the extent of the overdensity is comparable to the horizon: ro = 
At later times the perturbation becomes subhorizon. It is apparent that the form of the 
geodesic changes at a time f e = 1461, when the light enters the central homogeneous 
region of the overdensity. The endpoint of the geodesic is at f = at the present time 
fo = 1563. In the same figure we also plot the redshift of photons emitted at a certain 
time f em by sources located at the comoving coordinate f g {f em ) and observed at the 
present time fo at f = 0. The redshift of a photon emitted at the surface of the central 
homogeneous region is Zq = 0.012. 

In fig. 7 we depict the cosmological acceleration as a function of the redshift. The 
solid line is the acceleration parameter evaluated through eq. (5.14) for photons emitted at 
various points along the geodesic of fig. 6. The dashed line is the local radial acceleration 
given by eq. (4.1) at the point of emission. We observe strong deceleration for redshifts 
z < zq = 0.012. These correspond to emission points within the central homogeneous 
region that expands like a closed Universe. On the other hand, photons emanating from 
immediately outside the central region indicate strong acceleration. The acceleration 
remains positive up to redshifts z ~ 0.22. Asymptotically both acceleration parameters 
approach the value q = —0.5, typical of a flat, matter dominated Universe. 

In fig. 8 we repeat the calculation for a model with the same value e = 1.01, but with 
Hi = 4. This choice corresponds to a superhorizon initial perturbation, as r = 4/iTj. 
The present time is f = 1567 and the redshift of the surface of the homogeneous region 
Zq = 0.067. The central region is again strongly decelerating. However, there is a range 
of redshifts extending up to 0.5 or 0.7, for which q or q r , respectively, are positive. 

The direct comparison of the luminosity curves with the data from supernova 
observations is not feasible within our model. The reason is the presence of the strong 
deceleration for low redshifts. We have seen that at times f c ~ 1.55 -3 ' 2 = 1500 the 
central region stops expanding and reverses its motion. On the other hand, the absence 
of pressure in the Tolman-Bondi fluid does not permit a consistent description of the 
collapsing phase. Our model combines features characterizing the cosmological expansion 
at large scales and the growth of inhomogeneities at smaller scales. Its simplicity does not 
allow for quantitative accuracy in its predictions. However, our results demostrate the 
presence of a mechanism that could link the cosmological expansion with the appearance 
of large overdensities. 
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6. Summary and conclusions 

The notion of accelerating expansion has subtleties in inhomogeneous cosmology. The 
usual definition of the acceleration parameter in terms of eqs. (2.7)-(2.9) leads to eq. 
(2.10) for a pressureless and irrotational fluid. As a result the expansion is expected to be 
always decelerating if the energy density remains positive. It is not clear, however, if the 
definition of eq. (2.7) is the most appropriate for the characterization of the expansion. 
The Hubble parameter H in eq. (2.7) accounts for the volume increase at a given point 
during the expansion. In the presence of inhomogeneities or anisotropies H results from 
a certain averaging over the various directions. For the Tolman-Bondi metric this is 
obvious from eq. (2.14). On the other hand, the observation of light from a certain source 
in the sky provides information on the expansion along the direction to that source. 
Observing sources in various directions leads to a certain averaging. In inhomogeneous 
or anisotropic cosmologies it is not obvious that this is the same averaging that gives the 
volume expansion as quantified in the Hubble parameter H . 

If there are preferred directions in the underlying geometry one can define several 
expansion parameters. For one preferred direction, such definitions are presented in 
eqs. (2.12), (2.13), where the explicit expressions are also given for the case of the 
Tolman-Bondi metric. This metric describes the evolution of a spherically symmetric 
configuration. For an observer located away from the center of the configuration, the 
local Hubble parameter H and the acceleration parameter q are functions of the angle 
between the radial direction and that of observation. For signals originating near the 
observer, the parameters H and q can be evaluated by expanding locally the luminosity 
distance of an incoming signal in powers of the redshift. They are given by eqs. (2.16) 
and (2.18)-(2.22), respectively. In our study we considered the acceleration parameters 
q r , qe along the radial direction and perpendicularly to it. They are given by eqs. (2.23), 
(2.24), respectively. 

In the first part of this work (sections 3 and 4) we examined the effect of a mass located 
at the center of a spherically symmetric configuration on the surrounding cosmological 
fluid. Our cosmological solution is similar in spirit to the model of spherical collapse. 
We considered a fluid that is initially homogeneous and subject to uniform expansion. 
We assumed that at a certain time a spherical inhomogeneity appears at some point 
in space. We approximated the inhomogeneity as a point source located at the center 
of the configuration. Its presence modifies the local gravitational field and distorts the 
expansion of the surrounding fluid. For an observer located away from the center of such 
an inhomogeneity, the perceived local evolution is as follows: 

a) A central overdensity leads to acceleration along the radial direction, and deceleration 
perpendicularly to it. 

b) A central underdensity leads to deceleration along and perpendicularly to the radial 
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direction. 

c) A negative mass at the center leads to deceleration along the radial direction, and 
acceleration perpendicularly to it. 

In all cases the expansion becomes typical of a homogeneous, dust-dominated Universe 
far away from the inhomogeneity. 

Our solution demonstrates the link between gravitational collapse near overdensities 
and radial acceleration. Previous studies of expansion within the Tolman-Bondi model [1 
[12] have missed this point because they focused on the acceleration parameter defined 
through the volume expansion according to (2.7)-(2.9). This corresponds to a certain 
average of the expansion rates in various directions, as is obvious from eqs. (2.12)-(2.14). 
The acceleration defined in this way remains always negative, as expected from eq. (2.10). 
Also, in previous works the Tolman-Bondi metric has been studied in a gauge in which 
the initial energy density is constant in space, while the effect of the inhomogeneity is 
introduced through a function that determines the local Big Bang time. This obscures 
the intuition on the effect of large concentrations of mass on the expansion. 

It is apparent from our discussion that, if the observations lead to an effective 
averaging of the expansion rate over various directions similar to eq. (2.14), deceleration 
should be expected. If, however, the averaging is modified the expansion could be 
perceived as accelerating. This would happen in the case of an overdensity if the 
observations are more sensitive to the expansion along the radial direction. This is the 
case for an observer located at the center of an overdensity, who receives light signals only 
in the radial direction. In the second part part of this work (section 5) we studied the 
luminosity function for such an observer. The overdensity cannot be approximated by a 
point source any more, as the location of the singularity in the metric would coincide with 
that of the observer. For this reason, we assumed a continuous distribution of matter, 
that is homogeneous in the central region, falls off at larger distances, and asymptotically 
becomes homogeneous with an energy density smaller than the central one. We were 
mainly interested in the light signals that originate in the regions outside the central 
overdensity. 

The form of the luminosity distance as a function of the redshift in the above scenario 
provides another measure of the cosmological acceleration through eq. (5.14). For a range 
of redshifts corresponding to photons emanating from immediately outside the central 
overdense region, we observe strong acceleration. The acceleration remains positive up to 
a certain redshift, then turns negative, and asymptotically approaches the value q = —0.5, 
typical of a flat, matter dominated Universe. For small redshifts z ~ we observe strong 
deceleration. This range corresponds to photons emitted within the central homogeneous 
region, that evolves like a closed Universe and eventually collapses. We expect this feature 
to be modified in more realistic collapse models. These results are depicted in figs. 7, 8, 
and demonstrate the possible link between the growth of perturbations and the perceived 
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acceleration of the expansion. 

There are several points in the above scenario that require further work. One is the 
location of the observer, which is constrained by the isotropy of the cosmic microwave 
background. Exact isotropy would require the observer to be located exactly at the center 
of a spherically symmetric perturbation. Small deviations from complete isotropy, such 
as the difference of power between the northern and southern hemisphere or possible 
allignments of the low multipoles [21], could possibly allow for an off-center location. An 
alternative possibility is that the isotropy appears at a sufficiently large length scale in 
a configuration with several centers of symmetry. It is an interesting question whether 
the radial acceleration can be the dominant contribution to the averaged expansion for a 
random location of the observer. 

The form of the expansion within a collapsing overdensity is another point that has 
to be addressed in future work. In our model the cosmological fluid within the overdensity 
is approximated as pressureless. As a result, the expansion for small redshifts is strongly 
decelerating. This makes the comparison with the supernova data problematic. A more 
sophisticated model is required in order to describe with quantitative accuracy both the 
cosmological expansion and the collapse of an overdense region. 

Another issue concerns the existence of a characteristic scale for the inhomogeneity. 
In our simple model, it is apparent from figs. 7, 8 that the presence of acceleration at 
redshifts z ~ 0.5 forces the redshift of the surface of the overdensity to values zq ~ 0.05. 
On the other hand, it is remarkable that perturbations at the level of 10% in the present 
average density are capable of inducing acceleration in certain directions (see fig. 3). It 
is possible for such perturbations to be very extensive (of the order of 100 Mpc or even 
larger) without being in conflict with observations. 

As a final remark we mention that a series of studies [22] has considered the possibility 
that we are located within an underdense region of the Universe, while we receive light 
signals that originate outside this region. The luminosity-redshift relation in this case 
could be similar to the one in an accelerating homogeneous Universe. We have verified 
this conclusion within our model. The expansion is decelerating for all redshifts, but a 
significant reduction in the value of the Hubble parameter beyond a certain critical value 
z cr ~ 0.1 can provide consistency with the supernova data. 
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